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ABSTRACT 



We present a spectroscopic and photometric determination of the distance to 
the young Galactic open cluster Westerlund 2 using WFPC2 imaging from the 
Hubble Space Telescope and ground-based optical spectroscopy. HST imaging in 
the F336W, F439W, F555W, and F814W filters resolved many sources previously 
undetected in ground-based observations and yielded photometry for 1136 stars. 
We identified fifteen new 0-type stars, along with two probable binary systems, 
including MSP 188 (03 + 05.5). We fit reddened SEDs based on the Padova 
isochrones to the photometric data to determine individual reddening parameters 
Rv and Ay for 0-type stars in Wd2. We find average values (Ry) = 3.77 ± 0.09 
and {Ay) = 6.51 ± 0.38 mag, which result in a smaller distance than most other 
spectroscopic and photometric studies. After a statistical distance correction 
accounting for close unresolved binaries (factor of 1.08), our spectroscopic and 
photometric data on 29 0-type stars yield that Westerlund 2 has a distance (d) = 
4.16 ± 0.07 (random) -1-0.26 (systematic) kpc. The cluster's age remains poorly 
constrained, with an upper limit of 3 Myr. Finally, we report evidence of a faint 
mid-IR PAH ring surrounding the well-known binary candidate MSP 18, which 
appears to lie at the center of a secondary stellar grouping within Westerlund 2. 

Subject headings: open clusters: individual: Westerlund 2 



1. Introduction 



The open cluster Westerlund 2 (hereafter Wd2, Westerlund 1961), at 1 - 3 Myr, is one 
of the youngest massive stellar clusters known, similar in many respects to better-known 



^Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are associated with program # 8zt. 

^This paper includes data gathered with the 6.5 meter Magellan Telescopes located at Las Campanas 
Observatory, Chile. 

■^Based on observations obtained at the Southern Astrophysical Research (SOAR) telescope, which is a 
joint project of the Ministerio da Ciencia, Tecnologia, e Inovagao (MCTI) da Repiiblica Federativa do Brasil, 
the U.S. National Optical Astronomy Observatory (NOAO), the University of North Carolina at Chapel Hill 
(UNC), and Michigan State University (MSU). 
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clusters such as the Arches (Figer 2005) 



(Stolte et aL 2004 2006) 



Quintuplet ( |Figer et al.|[l999| , and NGC 3603 

(lO'^ 24" 



Located in the Carina arm at (a, 5) = (10'^ 24"* 01!1, — 
57° 45' 32"), (l, b) = (284?3, -0?3), it powers the surrounding giant H II region RCW 49 
(Rodgers et al. 1960). Figure [l] shows a three-color image of Wd2 and RCW 49 as seen in 
the Spitzer Space Telescope mid-IR Galactic Legacy MidPlane Survey Extraordinaire mosaic 
images (GLIMPSE, [Bi^n amin et al. I 20031). Blue represents the Infrared Array Camera 
(IRAC, Fazio et al. 2004) 4.5 /im band which highlights the reddened {Ay ~ 6) stellar 
cluster, while green (5.8 /xm) and red (8.0 fim) highlight the photodissociation region (PDR) 
where PAHs excited by soft UV photons re-radiate stellar energy at mid-IR wavelengths. The 
seminal imaging study on Wd2 by Moffat et al. ( 1991[ ) (hereafter MSP91) suggested that the 



cluster may contain more than 80 0-type stars. Churchwell et al. (2004) and Whitney et al 



(2004) used Spitzer mid-IR images to discover a complex network of dust pillars, filaments. 



and young stellar objects (YSOs) within RCW 49, suggesting that star formation may be 
ongoing, or even triggered in the surrounding clouds by winds and ionizing radiation from 
the central cluster. 



Dame (2007) and Furukawa et al. (2009) used millimeter- wave CO spectroscopy to 
conclude that the mass of the molecular cloud associated with Wd2 is 7.5 x 10^ Mq and 
1.7 ± 0.8 X 10^ Mq, respectively, sufficient to produce the massive stellar cluster Wd2. The 



only study that has tried to estimate Wd2's stellar mass is Ascenso et al. (2007). They found 



that the cluster's mass function is consistent with a Salpeter (1955) power law with an index 



of —1.20 ± 0.16. They estimated a total stellar mass of 7000 Mq, assuming a distance of 2.8 
kpc, but the mass may be larger if greater distances are adopted. If this is so, then Wd2 is 
among the most massive clusters in the Galaxy, similar to the Quintuplet cluster (6.3 x 10^ 



Mq, Figer et al. 1999), Arches (2 x 10^ Mq, Figer 2005), or NGC 3603 (7 x 10^ Mq, Stolte 



et al. 2004, 2006). 



Wd2 and RCW 49 have attracted considerable attention because of ongoing star for- 



mation (Whitney et al 2004; Churchwell et al. 2004), and the possible association with 



two very-high-energy (VHE) 7-ray source HESS J1023-575 and HESS J1023-5746, de- 



tected in the High Energy Stereoscopic System (HESS, Aharonian et al. 2007 Reimer et 



aL||2008 H.E.S.S. Collaboration et aL||2011 ), whose positions are coincident with the pulsars 



PSR J1023-5746, and PSR J1028-5819, respectively (Saz Parkinson et al. 2010 Ackermann 



et al. 2011; Abdo et al. 2009). This discovery has suggested possible new mechanisms for 



VHE emissions, such as 7-ray production in the colliding wind zone of WR 20a, collective 
effects of stellar winds from the Wd2 cluster, shocks and magneto-hydrodynamic (MHD) 
turbulent motion inside the hot H II bubble, or shocks driven by supernova explosions into 
the interstellar medium (ISM). This cluster is rich in X-ray sources, the most luminous ones 



being the very massive eclipsing binary Wolf-Rayet stars WR 20a and WR 20b (Tsujimoto 



et al. 


2007 


Naze et al. 


2008 



stars, implying the possibility of many close binaries in this cluster (Naze et al. 2008) 



Since the discovery of Wd2 there has been much controversy over the distance to this 
cluster. MSP91 used UBV charge-coupled device photometry and spectroscopy to obtain 



a distance of 7.9;*^J;o kpc. Piatti et al.| (1998) obtained spatially integrated optical spectra 



of the cluster and reanalyzed a subset of cluster O stars from MSP91 to infer a revised 
distance of d = 5.7 ± 0.3 kpc. Carraro & Munari (2004) analyzed new UBVRI photometry 
for the cluster, finding d = 6.4 ± 0.4 kpc and an age ^ 2 Myr, based on isochrone fitting. 



Using GLIMPSE photometry and spectra for a single 04 star, MSP 18, Uzpen et al. (2005) 
found a distance of ci = 3.23to;53 kpc. Analyzing optical spectra, light curves, and published 
photometry of the massive eclipsing binary system WR 20a ( WN6ha-|- WnNha) , along with 



stellar atmosphere models, Rauw et al. (2005) concluded that this system is a member of 
Wd2, giving it d = 7.9 ± 0.6 kpc. The apparent association between Wd2 and nearby 
molecular and atomic clouds led Dame (2007) to derive a kinematic distance of d = 6.0 ± 1.0 



kpc. Tsujimoto et al. (|2007) used Chandra X-Ray Observatory and near-IR imaging to 

2 Mq T-Tauri stars that serve as X-ray "standard candles" ; they 
2-5 kpc. Obtaining new photometry and new 



identify a population of 
obtained a loose distance constraint of d 



spectral classification of 12 0-type stars from MSP91 and the binary WR 20a, Rauw et al. 



(2007) (hereafter R07) found d = 8.0 ± 1.4 kpc. Comparing deep near-IR photometry in the 



(2007) found d 



J, H, and Kg bands with main-sequence and pre-main-sequence isochrones, Ascenso et al. 

2.8 kpc and an age of 2.0 ± 0.3 Myr. In one of the most recent works, 
2 — 1) observations to infer a kinematic distance of 

2 and ~8 kpc. 



Furukawa et al. (2009) used CO(J 



5.4lj^'4 kpc. In summary, Wd2 has been placed at distances between 



In view of this range of distances for Wd2, we have undertaken a program of optical 
photometry with the Hubble Space Telescope (HST) and spectral classification of additional 
MSP91 sources with the Magellan and the Southern Astrophysical Research (SOAR) tele- 
scopes in an attempt to resolve this discrepancy. In ^we describe the reduction of the HST 
WFPC2 imaging, photometry comparison with previous works, astrometry, and long-slit 
spectral data. In ^we analyze the color-magnitude diagrams (CMD), describe our spectral 
classification method, and individually comment on each observed source. In ^we analyze 
the HST photometry in conjunction with model isochrones and stellar spectroscopy to infer 
the reddening, distance, and age of the cluster. 
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2. WFPC2 Imaging, Photometry, & Optical Spectroscopy 

2.1. Observations 

2.1.1. HST 

The cluster Wd2 was observed using the Wide-Field Planetary Camera 2 {WFPC2) 
on the HST. The observations were performed during Cycle 13, proposal ID 10276, PI H. 
Kobulnicky, and HST program 8zt. A two-orbit HST imaging program was performed on 
2005 May 19 & 20 using the standard three-point line dither pattern with 1/3 pixel spacing. 
On the first orbit the cluster was centered on the PCI chip, and a series of three exposures, 
dubbed the "long exposures", were obtained in each of the filters F336W (160 sec), F439W 
(40 sec), F555W (8 sec), and F814W (20 sec). On the second orbit the cluster was centered on 
the WF3 chip, and a series of three exposures, dubbed the "short exposures", were obtained 
in the same four filters with exposure times of 50 sec, 12 sec, 2 sec, and 2 sec, respectively. 
Table [T] summarizes the observations. The solid and dashed polygons in Figure [T] depict the 
locations of the "short" and "long" exposures, respectively. The exposure times achieve a 5a 
detection at \^ = 20 on the "short" exposure images and at = 21.5 on the "long" exposure 
images. The UBV magnitudes of MSP91 were used to derive the exposure times to avoid 
saturation and bleeding of bright stars, but at the same time go as deep as possible down 
the mass function (MF) of the cluster. 



2.1.2. Optical Spectroscopy 

Longslit spectra were obtained for 10 MSP91 sources with the 6.5 m Magellan I Walter 
Baade Telescope and the 4.1 m SOAR telescope, using the Inamori- Magellan Areal Camera 
& Spectrograph (IMACS) and the Coodman High Throughput Spectrograph (GHTS), respec- 
tively. Owing to the crowding in the Wd2 core, other sources were caught serendipitously 
on the slit and augmented the number of observed sources to 34. 

Spectra for six stars (#528, #640, #714, #738, #857, and #906 in our numeration 
system, introduced in § 2.2.1) were obtained on 2011 April 6/7 with the Magellan Telescope. 
Each source was exposed for 300 sec along with a 5 sec exposure of NeHeAr calibration 
lamps with a 0'.'9 wide slit and a grating of 300 1 mm~^ which yields a reciprocal disper- 
sion of 1.34 A pixel"^. The typical seeing was 1" FWHM. The observations used 1x1 pixel 
readout modes yielding a scale of 0'.'2 pixel"^. The Magellan observations provided a wave- 
length coverage 4640 A < A < 10,113 A. Ultimately, only the spectral range from 4640 A - 
6581 A had a sufficient signal-to-noise ratio (S/N) to be usable. 
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The sources #137, #549, #1004, and #1031 were observed on 2011 June 9 & 10 
with SOAR. Typical seeing was I'.'S FWHM at SOAR on June 9/10, and T2 FWHM on 
June 10/11, as measured from the cross-dispersion width of the stellar profile. The charge- 
coupled device readout mode employed 2x8 pixel binning yielding a cross-dispersion scale 
of 1'.'2 pixel"^. Hence, the spectra are undersampled in the cross-dispersion direction. Three 
exposures of 60 sec each were obtained on each source and calibrated with 15 sec FeAr 
lamps. The 1'.'68 wide slit and the KOSI 600 1 mm~^ grating yielded a reciprocal dispersion 
of 1.3 A pixel" ^. The wavelength coverage includes the range 4600 A < A < 7300 A. 

We obtained spectra of stars #714 and #738 with the Magellan Inamori Kyocera Echelle 
(MIKE) spectrograph on the Magellan II Clay telescope on 2011 November 14 using a 1" x 5" 
slit in ~1'.'2 seeing. Exposure times were 2 x 600 sec for #714 and 1 x 900 sec for #738. 
The usable wavelength range covered 3900 A - 5050 A in 27 orders at a nominal spectral 
resolution of 0.09 A FWHM at 5000 A. 



2.2. Reduction 



2.2.1. HST Photometry 



We performed photometry on the long and short exposures in each filter using the PSF- 



fitting stellar photometry package HSTphot (Dolphin 2000) that was specially designed for 



the under-sampled PSF of the WFPC2. HSTphot provides both instrumental and standard 



UBVI magnitudes, using the standard reduction recipe and transformations of Holtzman 



et al. (2006) as updated by Dolphin (2009). This procedure includes masking out the bad 



pixels, generating sky images, masking cosmic rays, and cleaning hot pixels. HSTphot per- 
forms a two-pass search for stars in each WFPC2 chip to locate brightness peaks, performs 
PSF-fitting photometry in an iterative manner, and produces a "refined" photometric mea- 
surement and position for each star with an astrometric precision of 0.03 pixels. The code 
matches positions of stars detected in more than one filter to provide source-matched multi- 
color photometry. The initial run of HSTphot yielded 2251 and 577 candidate detections in 
at least one filter for the "long" and "short" images, respectively. 

It is a known characteristic of HST WFPC2 images that spurious sources arising from, 
for example, diffraction spikes near bright stars, are included in this list. In an attempt 
to remove false sources we applied a series of additional criteria to the initial source list. 
HSTPhot reports a series of "global" parameters for each source, based on the combined 
detections in one or more exposures and one or more filters. These include the source's 
roundness, sharpness, x (as a goodness of fit), object type (point source or extended), and 
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S/N. We determined a threshold roundness parameter, used to discriminate between point- 
hke and elongated objects, by fitting a one-sided Lorentzian profile to the distribution of 
roundness values and retaining only sources having roundness parameter less than la above 
zero, corresponding to roundness < 0.37 and < 0.19 for the "long" and "short" exposures, 
respectively. This retains 64% and 79% of the sources for the "long" / "short" exposures, 
respectively. The rejected stars are invariably those with the lowest S/N below ~10. Next, 
as per the HSTphot manual suggestion, we retained only sources having |sharpness| < 0.5. 
This criterion removes only a handful of stars that have not already been rejected by the 
roundness cut. Finally, we retained only sources with global S/N > 10 in order to ensure 
that the final catalog contains only highly reliable sources. Applying these criteria, the 
final list contains 967 and 360 stars in at least one filter for the "long" and "short" images, 
respectively: 194/171 in F336W, 227/185 in F439W, 348/286 in F555W, and 918/343 in 
F814W. We performed a visual inspection of the images to confirm that these selection 
criteria retained real stars while rejecting obvious spurious sources such as diffraction spikes. 

Figure [2] shows the photometric uncertainty as a function of instrumental magnitude 
for the short exposures (dots; colored blue in the electronic edition) and long exposures 
(asterisks; colored red in the electronic edition) in each filter. Typical uncertainties are 
< 0.01 mag for stars < 15 mag in both the long and the short exposures for all filters. 
Uncertainties rise to 0.2 mag at the faint end. This occurs for stars with F555W ^ 20 in the 
short exposures and F555W^ 21.5 in the long exposures. Stars with uncertainties exceeding 
0.1 mag appear in these plots because, on a per-filter basis, the S/N may be greater than the 
threshold global S/N of 10:1. The sudden drop in source density in the F814W panel (lower 
right) above uncertainties of ~0.1 mag arises because of the global S/N criterion imposed. 
Since F814W is least affected by reddening, most sources are detected more strongly in this 
band and, therefore, the F814W photon statistics drive the global S/N values. 

For the 190 sources detected on both the "short" and "long" exposure images, we 
compared the photometry in each HST filter to assess any systematic differences between the 
two orbits observations. Figure [3] plots the differences between the short- and long-exposure 
instrumental magnitudes for each HST filter as a function of long exposure magnitude. 
This figure shows that there are no systematic differences between exposures, except at 
the faintest levels where the asymmetric nature of Poisson noise produces a systematically 
brighter magnitude measured in the short exposures relative to the long exposures. However, 
the rms deviation from zero is much larger than the typical photometric uncertainty. The 
reduced exceeds 8 in all cases, indicating that there are additional sources of photometric 
uncertainty that are not included in the errors calculated by HSTphot. This is not an effect 
of chip position or crowding, but rather is a known limitation of existing HST photometry 



packages when dealing with very short exposures (Dolphin 2000). Therefore, we adopt 
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photometry from the long-exposure images, when available. For stars appearing on only the 
short-exposure images (138 of the 1136 stars), we adopt a larger photometric uncertainty 
that is the quadrature sum of the uncertainty estimated by HSTphot and an uncertainty 
determined empirically from the dispersion in Figure |3] using 1-magnitude wide intervals. 
This procedure yields xl values of less than two in all of the filters. Stars detected on the 
"short" or "long" exposure images are identified as such in the last column in Table [2j 

Coordinates for each star were obtained using the task METRIC in the Space Telescope 
Science Data Analysis System (STSDAS) package. For a subset of 120 bright isolated stars 
we identified counterparts in the 2MASS Point Source Catalog (PSC) and used the mean 
differences between positions to correct the headers of the HST images. The position header 
parameters obtained from the guide stars were corrected by (Aa, A6) = (0'.'27182, — 0'.'57479) 
for the "short" exposure images, while the "long" exposure images required a correction of 
(Aa, A5) = (— 0'.'01137, — 0'.'52219). After astrometric correction, the rms deviation between 
positions measured on the HST images and the 2MASS coordinates is O'.'l in R.A. and Dec. 
This is consistent with the limitations imposed by the WFPC2 pixel size coupled with the 
precision of the PSC astrometry. 

Table |2] lists the final band- and exposure-merged photometry in each filter, totaling 
1136 unique source IDs ordered by right ascension. A few lines appear in the printed journal 
as a sample of the table format. The entire content is available as a machine-readable table 
in the electronic edition. Column 1 gives the ID number from this work. Column 2 provides 
a cross-identification with the nomenclature of MSP91. In many instances the angular 
resolution of HST reveals several sources within a ~ 2" radius of the nominal MSP91 star. 
These are denoted by the MSP91 (MSP) designation with an alphanumeric suffix "a" for 
the brightest component, "b" for the next brightest component, as measured in the F555W 
band, etc. Columns 3 and 4 give the (J2000) Right Ascension and Declination in degrees. 
Columns 5 through 12 list the instrumental magnitudes and uncertainties in the HST filter 
system. Columns 13 through 16 list the transformed magnitudes in the standard UBVI. 
Column 17 lists the global S/N value. Column 18 lists whether the photometry is from the 
"long" exposure image (1) or the "short" exposure image (s). 

We caution that the transformations from the HST to Johnson system used by HSTphot 
are based on stars having a limited range of color, and that many of the Wd2 stars fall at the 
extreme red end of those calibrations. These reddest stars are late-type, low-gravity giants 
rather than heavily reddened early type stars with high surface gravity that typify the Wd2 



membership. Holtzman et al. (1995) further caution that the transformation are known 
to depend on metaUicty and surface gravity. We conducted our own synthetic photometry 
on reddened 0-star spectra to assess the extent to which the HSTphot transformations are 
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appropriate to our targets {{B — ~ 1.4; {U — B) 0.3.) We found that, for the B and V 
bands, the transformations are appropriate (within 0.02 mag) for reddened O stars. In the 
U and / bands, systematic uncertainties exceeding 0.05 mag may be present, as noted by 



Holtzman et ah (1995). Hence, where high levels of photometric accuracy are required, the 



HST instrumental magnitudes are preferred over the UBVI magnitudes. We use the HST 
instrumental magnitudes hereafter in our analysis, but we provide the URVI magnitudes as 
a convenience and compare them to the MSP91 values in the next subsection. 

Transformations from F336W to the standard [/-band are particularly problematic. 
One complication is the red leak in the F336 filter, which becomes more significant for 
redder stars. However, for the great majority of Wd2 stars in the range, —0.5 < (F336W - 
F439W) < 0.5, or {U-B) ~ 0.3 (MSP91), this only amounts to ~2-3%, as shown by Figure 3 



of Holtzman et al. ( 1995 ) and confirmed by our own synthetic photometry of reddened early 



type stars. Accordingly, we correct all of the reported U and F336W magnitudes by +0.02 
mag, acknowledging that this is only a zeroth-order correction appropriate to the mean 
Wd2 stars. More problematic is that the standard Johnson U band straddles the Balmer 



discontinuity while the F336W filter lies short ward of it. Holtzman et al. (1995) provides 



a discussion of difficulties in transforming the F336W magnitudes into standard [/-band 
magnitudes, notes that the F336W zero points and [/-band transformations are based on 
stars having {U — B) < —0.1 and {V — I) < 1.0, a regime that excludes the reddened O 
stars of Wd2. Accordingly, we urge caution in any interpretation of the transformed [/-band 
values. Nevertheless, our U-band photometry agrees well with that of MSP91. 

Figure |4] displays the positions of the 1136 detected stars in equatorial coordinates. The 
concentration of cluster stars located on the PCI chip of the "long" exposure images and 
in the WF3 chip of the "short" exposure images, is the dominant feature of the cluster, a 
secondary concentration of stars about 45" to the north is also apparent. Equally striking 
is the near-total absence of stars over a ~25" diameter area to the northeast of the Wd2 
core. This corresponds to a minimum in the mid-IR [5.8] and [8.0] emission in Figure [l| 
suggesting the presence of an infrared dark cloud (IRDC) having extremely high extinction. 



Our examination of the CO molecular maps of Furukawa et al. (2009) reveals no obvious 



molecular clouds at this location, although the angular resolution of those data is l!5, several 



times the size of the region of interest. The JHKg infrared images of Ascenso et al. (2007) 



show no obvious dearth of stars at this location, lending credence to the hypothesis of a 
localized region of high extinction that affects the optical wavelength data most significantly. 

Figures [5] - [7] show three-color WFPC2 images of Wd2 with F439W in blue, F555W 
in green, and F814W in red. The star ID numbers on the figures denote stars of interest 
to this study according to the numeration of Table [2] The color scale reveals that the core 
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cluster stars are red compared to other bright foreground field stars that appear white, 
such as MSP 158 (#436), confirmed as a foreground object by R07. The excellent angular 
resolution reveals that many of the bright core stars identified by MSP91 have close neighbors 
that are blended or unresolved at ground-based resolutions. 



2.2.2. Photometry Comparison 

We compared the photometry obtained with HSTphot to that presented by MSP91 and 
R07 for the 75 and 11 stars in common with those works, respectively]^ The catalog provided 
by MSP91 does not contain positions, so it was necessary to do a visual match between our 
images and their finding chart. Since the MSP91 observations were performed with ground- 
based telescopes that do not possess the resolution of the HST, some of the MSP stars are 
actually blends of several stars that are now resolved in our data. For such cases we choose 
the brightest star in the HST image to be the one that was originally identified by MSP91. 
In 63 of the 75 cases, either the MSP star is isolated and unambiguous, or we can identify 
the MSP source with the brightest of several close sources in our photometry. In 8 cases, 
we found two or more stars of roughly equal brightness within a radius of 2", so the cross 
identifications with MSP91 are ambiguous. 

Figures [8] and |9] show differences between our photometry and the works of MSP91 and 
R07 (Am = rriHST — fnthem) versus HST magnitudes transformed to the Johnson system. 
Numbered points indicate the nomenclature of MSP91 in both Figures. Where significant 
photometric differences exist, they are seen in both the V and B bands. The median dif- 
ferences compared with MSP91 are 0.06, 0.19 and 0.15 mag, for the U, B and V bands, 
respectively. Compared with R07 the median differences are 0.22 and 0.12 mag for the 
B and V bands respectively. Both figures show that our measurements are systematically 
fainter than theirs in both the BV bands, particularly for the brightest stars which lie in 
the most crowded regions. Part of this effect can be explained as a result of ground-based 
photometry that inevitably includes close neighbors in the measurement of the central star. 
These neighbors are well-resolved and individually measured in the HST data. 

To further investigate this systematic magnitude offset we performed a series of tests on 
the HST images. First, we verified that the images are processed using the latest calibration 
pipeline routine, which in this case is the OPUS version 2009_2k with calwp2 version 2.5.5 



■^Unfortunately, WR 20a (MSP 240) is saturated on the F555W and F814W images so no photometry 
is available in those bands. 
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(Apr 17, 2009)[^ This represents the final and best cahbration of the WFPC2. Next, we 
performed aperture photometry using the Image Reduction and Analysis Facility (IRAFJ^ 
package and compared it with the PSF results of HSTphot. To this photometry we applied 
the corresponding filter's zeropoint and Charge Transfer Efficiency (CTE) corrections ac- 
cording to Dolphin| ( |2009 ) and subsequent updates from the WFPC2 Calibration and CTE 
Corrections webpag^ We found excellent agreement between the results of the aperture 
and PSF photometry. Furthermore, the HST zeropoints are well calibrated and their un- 



certainties range between 0.02 - 0.04 magnitudes, depending on the filter used (Heyer et al. 



2004). 



Another possible source of discrepancy with ground-based observations could be the 
color transformations applied to the HST fiight system filters. The color transformation of 



Holtzman et al. ( 1995 ) was made using unreddened standard stars and stars from the cluster 



u) Cen rather than reddened O stars like those that dominate the population of Wd2. We 
used the IRAF/STSDAS synphot package to create synthetic photometry of reddened O 



stars and compare it with the color transformations of Holtzman et al. (1995). We found 



the transformations for reddened O stars to be compatible with the color transformations 
based on late- type stars, to within a few percent. In most cases the color transformations are 
indeed small and the HST^s fiight system filters can be used as a proxy for standard filters. 



given the Holtzman et al. (1995) transformations. For example, the differences between the 
F555W and the standard l^-band are on the order of ~ 0.01 mag for a wide range of stellar 



colors (see Figure 4 in Holtzman et al. 1995). 



We also investigated the known long-short exposure difference (Casertano & Mutch 



ler 1998) and the low background level (Golimowski & Biretta 2010) problems that affect 



the WFPC2 camera. Both effects are well-calibrated by HSTphot^ and any residuals in 
the current CTE corrections are small (~0.15 mmag/row for images with low-to-moderate 
background signal of < 50 e~) and cannot account for the systematic difference between the 
ground-based and WFPC2 observations. 

Others have also noticed discrepancies between HST and ground-based observations. 



For example, Dotter et al. (2008) reported offsets in the V and / bands of 0.068 and 0.007 



^An older version of the pipeline, OPUS version 2008_5c with calwp2 version 2.5.3 (Sept 4, 2008) produces 
the same photometry when using HSTphot. 

^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

^http:/ / purcell.as.arizona.edu / wfpc2_calib / 
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mag, respectively. 



Turner (2012) also reports offsets with UBV photometry. From his 
0.2 magnitudes. At the moment we can only conclude 



Figure 2 we estimated an offset of 
that there is an offset of unknown origin between the values reported here and those of 
MSP91 and R07. 



2.2.3. Spectral Reductions 

Optical spectra from Magellan and SOAR were bias-subtracted, flat-fielded, wavelength- 
calibrated and normalized using the standard IRAF procedures. The rms of the wavelength 
calibration was 0.09 A for the Magellan data and 0.28 A for the SOAR data. The spectral 
resolution determined from the width of arc lamp exposures was 4.1 A FWHM, measured 
at 6172 A, and 7.6 A FWHM, measured at 6172 A, for the Magellan and SOAR spectra, 
respectively. Signal-to-noise ratios of the reduced spectra range between 15:1 and 91:1 pixel"^ 
at 6000 A for the Magellan data and 11:1 to 80:1 pixel"^ for the SOAR data. For purposes 
of spectral classification, we used a subsection of the observed spectral range from 4600 A - 
66OOA having good S/N and a number of temperature-sensitive He I and He II lines. The 
blue spectral range (< 4600 A) commonly used for spectral classification was either not 
observed (Magellan) or unusable (SOAR) because of low S/N. 

The MIKE spectra of #714 and #738 were reduced using fiat fields from the internal 
quartz lamps with a diffuser screen. Spectra were extracted for each star and wavelength 
calibrated to an RMS of 0.004 A using ThAr lamps. Multiple spectra and spectral orders 
were combined, discarding data of low S/N within 10 A of the edge of each order, and 
continuum normalized using a 7th order Legendre polynomial before being Doppler corrected 
to the Heliocentric velocity frame. S/N ratios were 23 and 11 per pixel for #714 and #738, 
respectively, at 4500 A with 0.04 A pixels, but we smoothed the spectra to ~1 A resolutions, 
yielding S/N ratios of 85 and 50, respectively. 



3. Data Analysis 
3.1. Color Magnitude Diagrams 



Figure [10] shows the CMD for Wd2 in the HSTs standard filters. The CMD is the 
F555W magnitude vs. the (F439W - F555W) color. Detections with MSP91 counterparts, 
81 stars, are identified by an asterisk (colored red in the electronic edition). Stars of relevance 
to this work are labeled using the numeration system of this work. The CMD clearly show 
shows two distinct populations. On the CMD the blue "arm" represents the field stars, while 
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the redder "arm" clearly shows the cluster's members. Two stars, MSP 91 and 158 (#832 
and 436, respectively), fall on the field "arm" and are likely foreground objects, consistent 
with the conclusions of ROT. 

We attempted to clean the CMD of field stars by using the 3D field star decontamination 



algorithm developed by Bonatto & Bica (2007) and then improved in Bica et al. (2008) and 



Bonatto & Bica (2009). Unfortunately, the limited WFPC2 field of view means that we lack 



off-source observations of field stars needed to make a proper field star decontamination. 
However, given the richness of the cluster core in the center of the small WFPC2 field, these 
data are dominated by cluster stars and, at the brightest magnitudes, suffer little field star 
contamination in the red arm. 



3.2. Spectral Classification 



The vast majority of the new spectra for Wd2 exhibit lines of neutral and/or ionized 
He, indicative of O or early-B type stars. Four sources exhibit Ha in emission. Figures 11 



and [12] show the normalized spectra for early-type stars labeled according to the IDs of this 
work and organized by decreasing value of temperature/spectral type from top to bottom. 



Figure 13 shows the spectra for stars having spectral type consistent with late-O/early-B, 



undetermined, or field stars. The main stellar absorption lines, interstellar lines (mainly 
Na I D lines AA5889,5895), and Diffuse Interstellar Bands (DIB) are labeled. We classified 
stars having new spectra using the ratio of the equivalent width (EW) of He II (A5411) to 



the EW of He I (A5876) using the diagram and analytic fit from Kobulnicky (2012): 
EI^(A5411) 



1.16208 X lO^^^T^ - 1.19205 x 10^'T^ + 4.22137 x lO^^T - 50.5093. (1 



7m2 



EVr(A5876) 



This calibration is based on model spectra from the Lanz & Hubeny (2003) Tlusty models 



and hot stars from the Jacoby et al. (1984) spectral atlas. Because we have not observed 



spectral standards using our spectral setup, this may introduce additional uncertainty into 
our procedure for placing our targets on the MK system. We used the robust curve fitting 



package MPFIT (Markwardt 2008) as implemented in IDL to measure the EW of the He II 



A5411 and He I A5876 lines and their uncertainties. To improve the robustness of the fit 
for weak lines, we constrained the Gaussian width to that measured for the stronger of the 
two lines. Table |3] list the stars having spectral types determined using the above procedure. 
Columns 1 and 2 are the stars' identification numbers according to this work and MPS91, 
where available. Columns 3 and 4 are the photometry, V mag and B — V color, for such 
stars. Columns 5 and 6 are the measured EW for He II and He I lines having detections 
above 2cr; negative values correspond to emission lines, and upper limits are indicated. 
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Column 7 is the ratio of the He II/He I EW and column 8 is the stellar effective temperature 
corresponding to the ratio as determined by the relation in equation [T| Column 9 is the 



spectral type determined from the temperature approximation to one subtype (Martins et 



al. 2005). Because our SOAR/IMACS data do not cover the blue portion of the optical 
spectrum typically used for classification, we are not able to constrain the luminosity class, 
thus we adopt a default of dwarf (V). This is likely to be correct in most cases, since our 
targets are not among the brightest cluster members, assigned V or III by R07. A constraint 
on the luminosity is obtained when analyzing the MIKE spectra for stars #714, and #738. 



The lack of N III AA4634, 4640, 4662 in emission in spectra of Figures [16} and [17] confirms 
that these stars are not evolved. Remarkably, there are no known supergiants in the cluster, 
consistent with it being very young. 

#137 (MSP 165) — We tentatively classify this star as 04 on the basis of the relatively 
large He II EW (0.97 ±0.09 A) and relatively small He I EW (0.30 ± 0.07 A) in the SOAR 
spectrum (S/N = 50). This ratio indicates a very hot star with ~43,000 K, very close to 



the maximum range calibrated in the Kobulnicky (2012) diagnostic diagram. However, the 



detection of He I precludes classification as a true 03 given the Walborn (1971) criteria 
requiring its absence. This star is among the brightest in the sample at = 15.5 and 
also one of the reddest with B — V = l.QQ. It does not have detectable Ha emission. Its 
location on the cluster outskirts near a region of diffuse IR emission (see Figures [l] & [6]) is 
consistent with higher dust extinction at a location that has not yet been fully evacuated by 
the powerful winds of the central cluster. 

#260 (MSP 201) — Located on the cluster outskirts not far from #137 (see Figure [6]), 
we classify as a late O or early B because of the lack of a He II line from our SOAR 
spectrum (S/N = 21). Ha appears in emission, along with [O III] AA4959,5007, despite our 
best attempts at background subtraction, indicating the presence of small-scale variations 
in nebular emission near the star. 

#463 (MSP 125)— This star, together with #483 at 1'.'5 to the E, constitute MSP 125. 
Our SOAR spectrum of this star (S/N = 11) has insufficient S/N for spectral classification. 
With {B — V) = 0.87, this star is unusually blue for a cluster member. 

#505 (MSP 196) — This particular spectrum from SOAR (S/N = 27) is dominated by 
light from #505 but may contain contributions from #495 and #502 located within about 
1". Hence, MSP 196 is a blend of at least three stars. The dominant component, #505, has 
{B — V) = 1.3 and spectral type 08.5, consistent with cluster membership. 

#528 (MSP 229) — This source and its companion #523, separated by ~0'.'4, constitute 
MSP 229. Our two spectra of this object, one from Magellan (S/N = 42) and one from 
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SOAR (S/N = 53), yield 08.5 and 08, respectively, so we adopt 08 as a weighted mean of 
the brighter component which dominates the spectrum. 

#547 (MSP 28) — This V = 17.4 source has [O III] emission and Ha emission. He II is 
a non-detection, and He I is marginally detected in our S/N = 22 SOAR spectrum, making 
this a late-0 or early-B star. The position of this star is interesting since it is located near 
the secondary cluster, centered on MSP 18, and it lies upon a ring of diffuse mid-IR emission 



which surrounds the secondary cluster. Figure [M] is a Spitzer 3-color image of the region 
surrounding MSP 18 showing this ~ 36" diameter ring-like structure, which is present in all 
IRAC bands but best seen at [8.0], consistent with PAH emission from a photo-dissociation 
region. We believe that this is the first recognition of this feature, which possibly demarcates 
the boundary of the H II region surrounding MSP 18, although MSP 18 does not lie at its 
geometric center. 

#548 (MSP 151) — Although this star is located ~13"NW of the cluster core in an area 
that may be affected by small-scale dust extinction. Observed with SOAR (S/N = 80), a 
spectral type of 04 was determined. This classification is at odds with MSP91 and R07 
classification of 06 - 7III. The He II is very strong with EW = 0.94 ± 0.04 A and He I is 
rather weak with EW = 0.31 ± 0.06 A. Nebular lines are well subtracted, so we do not have 
an explanation for the earlier spectral type compared to R07. 

#549 (MSP 44) — Star #549, along with the faint companion #538 constitute the 



source MSP 44. Indeed, Rauw et al.| (hereafter Rll, 2011) identify this system as a P = 
5.176 d eclipsing binary (AV^ = 0.6 mag) composed of an BIV -|- pre-main-sequence star, or 
possibly even a triple system. We obtained two spectra for this source, one as target (S/N = 
59), and the second being serendipitous (S/N = 28). Both observations, obtained with the 
SOAR telescope, yield consistent spectra where the He lines are very weak or non-detections. 
For this star we also detected emission at [O III] and Ha. The presence of emission lines 
may indicate a localized knot of nebular emission. The residual He I emission may affect 
the stellar EW measurement, rendering stellar classification unreliable, but consistent with a 
very late-0 or early-B star, in agreement with Rll. Rll reports nebular emission associated 
with MSP 18, MSP 223, and several other stars. Evidently, the nebular structure in the Wd2 
region is complex both spatially and, as indicated by Rll, in velocity as well. Like #547, 
which also exhibits nebular emission, #549 lies near the ring of mid-IR emission encircling 
MSP 18. Mid-IR PAH emission usually traces the UV-heated surfaces of molecular clouds. 
Such a location could indicate that the MSP 44 system is very young, consistent with the 
interpretation of Rll that the secondary star is a pre-main-sequence object. 

#556 — Located ~66" to the north of the cluster core, #556 {V = 17.1) seems to be 
surrounded by diffuse IR emission, very likely coming from the ring surrounding MSP 18. 
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He I is detected in our S/N = 23 SOAR spectrum with an EW of 0.53 ± 0.11 A, but He II 
is not detected, yielding only an upper limit on the spectral type of late O or early B. 

#55^ (MSP 157b) — Located at ~10" from the cluster core where blending becomes 
problematic, MSP 157 contains contributions from as many as seven stars, the brightest 
three of which are stars #597, #584, and #568 located within ~ 1'.'2 radius. Our Magellan 
spectrum (S/N = 58) and classification implies an 08. This is "slightly" later type than 
the 05.5 - 6.5V by R07, but their spectrum is likely dominated by the brighter component, 
MSP 157a (#597). 

#620 (MSP 96) — Located in a relatively isolated region at ~26" to the north of the 
main cluster, #620 along with #663, constitute MSP 96, and its photometry is consistent 
with the measurement from MSP91. Both He lines are non-detections in our S/N = 43 
SOAR spectrum. This limits our classification to a very late O or early B star. 

#64-0 (MSP 233) — This V = 16.2 source is relatively isolated and lies close to diffuse 
IR emission at ~ 14" to the SW of the main cluster. The clear detections of the He lines in 
our Magellan spectrum (S/N = 36) provide a classification of about 09.5. 

#664 (MSP 188) — This bright star has two faint (AV" = 4.3 mag) neighbors at ~1'.'5. 
The very strong He II and the weak He I in our S/N = 91 Magellan spectrum yield an EW 
ratio of the He lines of 8.99 ± 4.20, placing this star at the upper limit of our classification 
scheme; therefore, we type it as an 04 or earlier. This is in agreement with R07 who 
revised the MSP91 spectral type to 04. As portrayed in the velocity plot of the He II 



and He I lines in Figure [15] the spectrum shows that both He lines are double-peaked, with 
He I showing two approximately equal-depth absorption lines, and He II displaying different 
strength absorption lines. We measure the velocity difference between the peaks for He I to 
be Aw ~ 360.5 ±103. km s'^and for He II to be ~ 489.6 ±36.8 km s~^ The He lines were 
deblended to calculate individual EW and EW ratios, so as to determine the spectral types 
of what we infer to be a binary system. The two components have EW ratios of 5.70 ± 3.32 
and 2.04 ± 1.18, implying spectral types of at least 04 and 05.5, respectively. The binary 
nature of this system partially explains why #664 is so luminous compared to other very 
early members of Wd2. 

#704 (MSP 175)— Star #704, along with the faint and close neighbor #713 (lying 
~0'.'44 away, midway between #704 and #714), forms MSP 175. It shows very strong He II 
in our S/N = 47 Magellan spectrum, while He I is very weak (EW = 0.15 ± 0.05), implying 
a spectral type near 04. This is slightly earlier than the 05V - 06V type assigned by R07. 
R07 notes that the He lines in MSP 175 are broad making classification uncertain; He II 
lines appear asymmetric with a blue wing in Figure 3 of R07. Our Magellan spectrum shows 
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asymmetry with a red wing in both the He II and He I lines, suggesting a possible binary 
system. 

#7i^ — Star #714, located at ~3"to the W of the main cluster core, is in a region of 
heavy crowding. Stars #664 (MSP 188), and #704 (MSP 175) lie in a crowded region near 
the cluster core within I'.'S of several other stars, including #714 (no MSP identification) and 



#713. Figure 16 shows the MIKE spectrum of #714 (S/N = 91), smoothed to a resolution 
of about 1 A with key spectral features of hot stars labeled. This spectrum is consistent with 
the lower resolution IMACS spectrum in that He I A4471 is extremely weak (not detected 
here), while the He II lines are very strong, indicating an extremely hot star. By comparison 



with the 0-star atlases of Walborn & Fitzpatrick (1990) and Walborn et al. (2002) we assign 



this star an 03 spectral type and dwarf (V) luminosity class. Signatures of 02 spectral 
types, such as N V AA4604,4620 absorption or strong emission lines of N IV A4058 and N III 
AA4634, 4640, 4642 are not seen, although these latter N features are weakly present. The lack 
of He emission lines and lack of strong N emission lines further supports a dwarf luminosity 
classification. Nebular [O III] AA4959,5007 and H/3 lines also appear in this spectrum, but we 
consider these to originate from the diffuse ionized gas that pervades the RCW49 H II region; 
these emission lines remain after background subtraction, presumably because of small scale 
variations in the nebular emission along the slit. In our single spectrum we find no evidence 
of binarity. However, all of the principle spectral stellar features in #714 are redshifted by 
about 1 A (67 km s^^) compared to expected rest wavelengths, and compared to molecular 



gas associated with the cluster which spans the range -11 to 11 km s ^ LSR (Ohama et al. 



2010). This redshift may constitute evidence of being a single-lined spectroscopic binary or 



may indicate that #714 has a peculiar velocity sufficient to unbind it from the cluster. 



#138 (MSP 168)— Figure [T7| shows our single MIKE spectrum of #738 {V = 14.9) 
boxcar smoothed to about 1 A resolution (S/N = 48). He I A4471, with an equivalent 
width of 0.29 A, appears in this spectrum, making #738 later than 03. The EW ratio 



EW(A4686)/EW(A4471) is 1.86, similar to 05 - 06 stars from the Walborn fc Fitzpatrick 



(1990) atlas. The radial velocity, as judged from the mean of key spectral features, is 



consistent with zero km s ^ LSR, making it consistent with the molecular gas surrounding 



the cluster (Ohama et al. 2010). We adopt a classification of 05.5. 



#769 (MSP 219) — This star is located at ~ 9" to the S of the main cluster core in a 
relatively isolated area. He II is a marginal detection in the Magellan spectrum (S/N = 31) 
but He I is strong with EW = 0.85 ± 0.05. This suggests that #769 is 09.5 or later. 

#771 (MSP 167)— MSP 167 is a blend of the stars #804, #771, #765, #777, and 
#772. A broad He II line and a well-defined He I line produced an EW ratio of 1.00 ± 0.13 
implying that this star is 08. 



- 18 - 



#857 (MSP 444) — MSP 444 is one of those hard-to-identify stars owing to its position 
in the cluster core and its close proximity to #843 (MSP 203; r^Qf'.Q). With HST MSP 444 
was resolved into three components, #857 and #895, and #824. The Magellan spectrum 
(S/N = 78) of #857 possesses a strong and well-defined He II line while the He I is broader 
and weaker than He II, implying a hot star of type 04.5, consistent with 04 - 05 by Rll. 
With {B — V) = 1.4, its color and magnitude are about right to be an early type cluster 
member. 

#878 — Located at ~8"to the SW at the outskirts of the main cluster, the two Magellan 
spectra obtained for this faint star {V = 16.6) have insufficient S/N (S/N = 15) to perform 
a reliable classification. He II is not detected, and He I is weak, hinting at a late O or early 
B. 

#879 (MSP 235) — This star is in relative isolation at ~16" to the SE of the main 
cluster. The Magellan observations, with a S/N = 36, show strong He I but minimal He II, 
implying a spectral type of 09.5. 

#896 & #903 (MSP 183 subcomponents) — These stars are located on the outskirts 
of the main cluster core where heavy blending in ground-based observations is inevitable. 
Within a radius of ~1'.'5 there are seven stars that could contribute to the brightness of 
MSP 183, namely, #869, #896, #903, #826, #880, #887, and #856. These two stars, 
separated by ~ 1'.'5, became partially blended on the longslit observations with the Magellan 
Telescope (S/N = 43 for #896 & S/N = 39 for #903). The strong He lines in both stars 
resulted in an EW ratio of 0.83 ± 0.13 for #896 and 0.90 ± 0.14 for #903, implying an 08.5 
type for both. 

#1004 (MSP 32), #1012 (MSP 24), and #1026 (MSP 20)— This stellar trio, observed 
spectroscopically with SOAR, is located on the diffuse mid-IR ring that surrounds the sec- 
ondary cluster, at about ~ 12" to the W of MSP 18. Star #1004's spectrum (S/N = 53) 
possesses strong He lines, although He II is noisy and broad, implying a classification of 
09.5. The spectra of both #1012 (S/N = 36) and #1026 (S/N = 14) display [O III] and Ha 
in emission and have insufficient S/N for reliable classification. 

#1028— With V = 14.3, and {B - V) = 0.37, this star is clearly a field star in the 
vicinity of the cluster core. The serendipitous SOAR spectrum (S/N = 62) shows only strong 
Balmer lines with EW(Ha) = 7.67 A, making it a probable A star. 

#1031 — Greatly isolated from the cluster core at ~73" to the SE it lies projected against 
a bright mid-IR rim illuminated by the central cluster. The S/N = 32 SOAR spectrum is 
heavily contaminated by emission lines, and no spectral type determination is possible. 



Interestingly, as shown in Figure 13, #1031 shows Mg II A5173 A absorption lines. 
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4. Results 

4.1. De- Reddening 

We measured the extinction, Ay, and the ratio of selective to total extinction Ry = 
Av/E{B — V), for each star by individually fitting a spectral energy distribution (SED) 
based on the known spectral type/temperature of the star and the theoretical absolute 
magnitudes intrinsic colors of the Padova stellar evolution tracks and isochrones (CMD 
v.2.2)0 dGirardi et aLl|2000| [Marigo et all2008| [Girardi et all2002| ). Each SED was reddened 



by each of three reddening laws ( Cardelli et al. 1989 ; Fitzpatrick 2004 ; Fitzpatrick & Massa 



2007) (hereafter CCM89, F04, and FM07 respectively) and compared, after appropriate 

^and Ascenso et al. (2007) JHKg photometry to find a global 

and Ay- 



normalization, with the HS 



chi-squared minimum over all plausible values for Ry, and Ay. We note here that, while 
CCM89 provide an analytical approximation of the mean interstellar extinction with Ry and 
Ay as free parameters, FM04 tabulates 38 reddening parameterizations as a function of Ry, 
and FM07 tabulates parameterizations of 328 sightlines, creating 328 discrete reddening laws. 
As an additional method to estimate Ry we also used the method of [Fitzpatrick fc Massa 



(hereafter FM09 2009) which is based on an average relation between BVKg photometry 



in their Figure 7. The fitted methods are sightly dependent on the anchor point used to 
normalize the Padova magnitudes to the data. To eliminate this dependency the fitting was 
done using each of the seven photometric values, when available, as anchors points and the 
resultant values of Ry, Ay and Ax are weighted averages of all seven fits. Figure 18 shows 
examples of the best fits for stars #137, 528, 704, and 1039. Each subplot displays the HST 



photometry as squares and the Ascenso et al. (2007) IR photometry as circles. The lines are 



the best fit to the photometry after reddening the predicted magnitudes for the appropriate 
star's spectral type. For comparison purposes the best fit for each method was included: the 
solid lines is the CCM89, the dash line is the F04, and the dotted line is the FM07 reddening 
laws (each one is colored in blue, cyan and magenta in the electronic edition). 

Table |4] compiles the fit results and shows that best-fitting values range between 5.71 < 
Ay < 7.54 and 3.44 < Ry < 4.31. Column 1 lists the star's ID according to this work, 
column 2 lists the adopted spectral type, columns 3-6 are the best-fitting Ry using the 
methods of CCM89, F04, FM07, and FM09 respectively, and columns 7-10 are the the 
best-fitting Ay values. The results indicate that Ry is significantly higher than the mean 



'Available from the website: 



Holtzman et al. 



littp://stev. oapd.inaf.it/cgi-bin/cmd_2. 2 



( 1995 ) caution that reddening determinations should be made in the HST filter system 



rather than the transformed Landolt UBVI system. 
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interstellar value of Ry = 3.1, meaning that the obscuring dust produces more extinction for 
a given amount of reddening, often interpreted as a result of a population of larger "gray" 
dust grains. On the other hand FM07 and FM09 argue that there is no universal reddening 
law or extinction curves and that extinction cannot be parameterized in terms of a single 
variable such as Ry- The so-called "anomalous" dust would, in reality, be "normal" for the 
Wd2 region and may be different from any other dust in the Galaxy. 

Stars #664, 869, 889 are not included in Table |4] since these only have photometry in 
three bands (F336W, F439W, F555W) because of saturation or blending. We regard that the 
adopted fit procedure provides unreliable results and we adopt, for these stars, the average 
Rv and Ay from the 29 other stars. Excluding these three stars, the average values for the 
cluster are (Ry) = 3.77 ± 0.09 and (Ay) = 6.51 ± 0.38 mag where the uncertainties are 
the dispersions of the sample. We conclude that the mean Ry is substantially larger than 
mean interstellar values, and that all stars are consistent, within the uncertainties, of this 
mean. On the other hand, there appears to be a significant scatter in Ay, consistent with 
appreciable internal extinction variations. 

Table [5] compiles the corresponding Ax for the HST filters corresponding to the best 
model fit in Table |4| Column 1 lists the star's ID according to this work, column 2 lists 
the adopted spectral type, columns 3, 4 & 5 are the Apasew values for CCM89, F04 and 
F07, and column 6 contains the weighted average of columns 3, 4, & 5. The same pattern is 
followed by columns 5 thru 18 for filters F439W, F555W, and F814W. We use the values of 
Table [5]to deredden the listed stars in the HST filter system and correct for the heavy, patchy 



extinction toward Wd2 which produces considerable dispersion in the raw CMD (Figure 10). 



The (F336W - F439W) vs (F439W - F555W) color-color diagram (CCD) in Figure [19 
shows a portion of a reddened Zero- Age Main Sequence (ZAMS) Padova isochron^^ using the 
global A\ averages from Tables |4] & [s] (Ry = 3.77, and Ay = 6.51 with corresponding color 
excess of E(F336W - F439W) = 1.47 and E(F439W - F555W) = 1.81). The dashed line 
(magenta in the electronic edition) is part of the reddening vector connecting the unreddened 
ZAMS O stars and the reddened Wd2 stars. Asterisks (red in the electronic edition) represent 
Wd2 stars having spectral types that we used to obtain the reddening parameters. The 
reddened isochrone provides a good match to the cluster stars, given that these stars exhibit 
a dispersion that is parallel to the vector, indicating substantial differential reddening and 
consistency with the mean reddening parameters found for this sightline. Additionally, the 



10 



The Padova isochrones are based on the non-overshoot models of ATLAS9 by 



CasteUi et al. 



(19971, 



which in turn recommends using Martins et al. ( 2005 ) observational values for O stars to transform effective 
temperatures and surface gravities into spectral types. We adopt the Martins et al. (2005) observational 
calibration. 
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dotted curve (cyan in the electronic edition), is a reddening track that takes into account 
small systematic color trends due to the broadband nature of the HST filters and the large 
reddening towards Wd2. The curve was constructed by performing synthetic photometry 
with the software package synphot using a CCM89 reddening law. As pointed out by 



Schmidt-Kaler (1981) and Turner (1989), the effective wavelength of a filter changes as a 



function of spectral type and reddening, leading to curved reddening tracks. Along this curve 
we label points that correspond to particular extinction values Ay. According to Turner 
(1989) the difference between a linear and a curved reddening track is ~2% per E{B — V). 



Although this is a small factor we included it in Figure [19] for illustrative purposes. When 
compared to the random and systematic errors of the intrinsic and observed colors that 
are discussed previously and in the next section, this difference is small so we adopt the 
straight vector in our reddening analysis. From the figure it is evident that there exists 
a systemic offset between the linear vector and the curved track, but this is small when 
compared to the dispersion in the photometry. This offset can be explained by the method 
we used in the creation of the linear vector. This vector was created using the average Ax 
from the reddening laws of CCM89, FM04, and FM07. Such values do not necessary have 
to reproduce the curve created by synphot that only uses the CCM89 reddening law. 

Support for the "anomalous" value of Ry toward Wd2 can be found in the studies of 
other open clusters and their Cepheid stars in the vicinity of the Carina arm. The open 
cluster Ruprecht 91, CoUinder 236, and Shorlin 1; with projected distances of ~3?2, ~5?6, 
and ~6?4, respectively, away from Wd2, also have different values from the nominal Ry = 3.1; 
Ruprecht 91 with Ry = 3.82 ±0.13 ( [Turner et al"][2005| ), CoUinder 236 with Ry = 3.82 ±0.13 



(the Ry was an adopted value because of similar reddening law with Ruprecht 91; Turner 
et al.||2009[ ), and Shorlin 1 with Ry = 4.0 ±0.1 ( [Turner[[20l2| . Furthermore [Turnerj ( [20l2j ) 



made a compilation showing that various clusters in the Carina complex have values greater 
than 3.1 in general (see his Table 1). 



4.2. Distance Determination 

We adopt the Padova isochrones to determine the cluster's distance modulus (DM). 
Isochrones with ages between 1-4 Myr and solar metallicity {Z = 0.019) were chosen to 
match the expected properties of Wd2. The Padova database provides absolute magnitudes 
(and therefore intrinsic colors) for stars in the HST filter system which we use to compute 
a spectroscopic and photometric distance for each star, following the prescription used by 



Hanson (2003) in finding the distance to Cygnus 0B2. We adopt either our own spectral 



classifications or those of R07 and Rll (based on the Walborn & Fitzpatrick (1990) call- 
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bration) in order to calculate the appropriate ^F555W ^^"^ (F439W - F555W)o values and 
their associated uncertainties by interpolating the Padova data tables. For stars typed using 
our procedure described above, we estimate absolute magnitude and color uncertainties by 
propagating errors from our determination of the effective temperature. For stars typed 
by R07, we adopt temperature uncertainties given in that work. We deredden each star 
according to the values in Table [5] and then solve for the distance modulus of each. 

Table |6] displays the results of our distance analysis. The columns list (1) the star ID 
numbers according to this work, (2) the MSP91 identification numbers, (3) the star's spectral 
type, (4) the observed apparent F555W magnitude, (5) the observed (F439W - F555W) 
color, (6) the absolute F555W magnitude, (7) the intrinsic (F439W - F555W) color, (8) 
the extinction ^P555W, (9) the distance modulus, and (10) the distance, in kpc, obtained 
for each star. The weighted mean are {DM) = 12.98 ± 0.03, and (d) = 3.85 ± 0.06 kpc 
(the simple mean yields (d) = 4.20 ± 0.06 kpc). This weighted mean is taking into account 
uncertainties in the spectral/luminosity classification. We exclude from this calculation stars 
#664, 869, 889, which, as stated in the previous section, lack complete photometry. This 
distance is substantially smaller than the estimate that would result if we were to assume 
the canonical Ry = 3.1. If we were to assume such, along with (Ay) = 5.16 (calculated from 
the Ay values in R07), then the resulting distance modulus becomes (DM) = 14.246 ±0.038 
and {d) = 7.07 ± 0.12 kpc, about twice the distance than when using {Ry) = 3.77. 



Figure 20 shows the extinction- and distance-corrected Wd2 CMD in the HST filter 
system, illustrating the application of the results from Table [6| along with ZAMS, 1-4, and 
10 Myr isochrones. Each star was individually dereddened, extinction corrected, and then 
shifted by the weighted mean DM. For easier comparison, key spectral types are labeled 
at the right edge of the figure. Each star is labeled with its ID and spectral type. The 
dispersion around the ZAMS is minimal, and very few stars have an observed spectral type 
inconsistent with the spectral type expected along the ZAMS isochrone. Star #889 (MSP 18) 
is not shown in the figure because of its odd position in the CMD after adopting the average 
Ax (overcorrection in color). Since this star lies to the north of the main cluster, it is very 
likely that the average reddening determined for the main cluster is not applicable to this 
star. 

The presence of close, unresolved binaries in a stellar population produces brighter 
apparent magnitudes than if no binaries were present. Consequently, the derived distances 



are systematically low compared to the single-star case. Kiminki et al. (2012) modeled 
the light contribution from secondary components among massive binaries in Cygnus 0B2 
and found that they contribute 16% of the luminosity, on average. This translates into a 
systematic distance error of 8%, in the sense that the derived distances are too small unless 
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binaries are taken into account. We apply this statistical correction to find a distance of 
4.16 ±0.07 kpc. 

A possible source of systematic error in our distance measurement is the intrinsic colors 
adopted for massive stars. There is a difference in the intrinsic colors of stars between the 
Padova isochrones and Martins & Plez (2006) colors. The Martins & Plez (2006) {B — V)q 



colors are, on average, 0.04 mag redder than those that result from the Kurucz models or the 



photometric measurements of Johnson ( 1966 ) and Fitzgerald ( 1970 ). This difference implies a 



smaller color excess and, therefore, a smaller extinction and larger distance modulus. [Martins 
& Plez (2006) warn that a redder intrinsic color introduces a difference of AAy ~ 
3.1. For our mean value of Ry 



0.124 mag 

3.77, the change in extinction required is 



assuming Ry 

AAv = 0.15, equivalent to a systematic distance error of +0.26 kpc. We conclude that our 
measurement of the Wd2 distance is 4.16 ± 0.06 (random) +0.26 (systematic) kpc. 

The distance that we determined favors a small value, like some in the literature, but not 



as small as the 2.8 kpc suggested by Ascenso et al. (2007). Ascenso et al. (2007) measured 



a distance by fitting deep JHKg CMDs to pre-main-sequence isochrones after correcting for 
reddening, assuming a single value for Ay derived from the four brightest cluster stars. De- 
spite adopting an Ay smaller than our preferred value, they obtain a smaller distance, which 
we attribute to uncertainties in the pre-main-sequence models and the inherent ambiguities 
of fitting such models to the data. Furthermore, the distance calculation is based on a prede- 
termined age for Wd2, which has not yet been properly constrained. Our measurement, like 
that of Uzpen et al. (2005), falls in the middle of the d = 2 - 5 kpc constraints determined 



in the X-ray study of pre-main-sequence stars in Wd2 (Tsujimoto et al. 2007). Our result is 



also very close to the 4.2 kpc kinematic distance adopted by Churchwell et al. (2004) on the 



basis of 21-cm absorption arguments (Caswell & Haynes 1987). Other kinematic estimates 



from the CO radial velocity of the associated molecular cloud (Dame 2007 Furukawa et al. 



2009) are slightly larger, 5-6 kpc — still broadly consistent with our result but with larger 



uncertainties owing to the range of molecular cloud velocities along this sightline. However, 
our measurement is inconsistent with distances larger than 5 kpc. The main reason why 



our preferred distance is smaller than that found by Piatti et al. (1998); Carraro & Munari 



(2004), R05, R07, and Rll is that all of these works adopted a standard Ry = 3.1 which 



leads to a substantially larger distance. Our adopted distance of 4.16 kpc places Wd2 com- 
fortably within the Carina spiral arm and near the solar circle, about 40% further away than 



where it is pictured in Figure 3 of the review of massive stellar clusters by Portegies Zwart 



et al. (2010). 



Due to the discrepancy between our transformed UBVI photometry and that reported 
via ground-based observations, as attested in § |2.2.2[ we consider here what would happen 
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if we brighten our flight system photometry by 0.2 magnitudes to achieve better agreement 
with those works. When such is the case, we find that Wd2 has a shghtly smaller reddening, 
(Rv) = 3.60 ± 0.10 and (Ay) = 6.23 ± 0.39. This smaller Ay more than offsets the 0.2 mag 
photometric correction and yields (DM) = 13.26 ± 0.03, and (d) = 4.38 ih 0.07 kpc. After 
corrections for binary stars this becomes 4.73 ih 0.08 (random) +0.26 (systematic) kpc. In 
summary, the photometric offset compared to other works has only a small impact on the 
final distance determination. This alternative distance is still much smaller than the "far" 
distances in the literature, substantially larger than some of the "near" distances, and still 
broadly consistent with the kinematic distance estimates. 



4.3. Cluster Age 



In Figure [20] we used the values of Table [5] to deredden all the available photometry 
from Figure 10, and we compared the data with 1-4 Myr isochrones in order to constrain 
Wd2's age. From the inspection it becomes evident that isochrones with ages 1-3 Myr are 
degenerate, and it can be "safely" assumed the majority of the stars are dwarf class. The 4 
Myr isochrone shows that there should be a population of evolved giant stars that is absent 
in the data. At the most we can conclude that Wd2 has an age between 1 and 3 Myr. Other 
authors have determined the same age constraint on the basis of integrated spectra, CMD 



fitting, and the lack of evolved stars within the cluster (Piatti et al. 1998; Ascenso et al. 



2007). 



We attempt here a different kind of age constraint using the size of the cavity carved in 
the ISM by the winds of massive stars ( Churchwell et aL||2004 ) . Weaver et al. (1977) derives 
(Equation 51) the radius of an interstellar bubble as a function of time in Myr, ambient 
density no, and mechanical energy injection L^q 
expression to find 



in units of 10 erg s . We invert this 



27 



1 1/3 



no 



Myr. 



(2) 



From the angular size of the cavity in Figure [T] (see also Churchwell et al. 2004 ) we adopt 



a radius of 105" or 2.0 pc at the adopted distance. We assume no = 10 cm ^, typical 
of molecular clouds. Based on the wind speeds and (highly uncertain) mass loss rates for 



O stars from Mokiem et al. (2007) we estimate a time-averaged energy injection rate of 
L36 = 1.0 from the ensemble of massive stars. The mechanical wind luminosity is dominated 
by the few most massive stars, so the exact number and type of stars later than about 05 
is inconsequential. The resulting expansion timescale of the (still ill-defined) interstellar 
cavity is ~0.6 Myr, indicating extreme youth, consistent with the abundance of molecular 
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material just outside the cluster core. We caution that the Weaver et al. (1977) analysis was 
developed for interstellar bubbles in uniform media and densities more typical of the diffuse 
ISM. The presence of a radio "blister" structure on the west side of the Wd2 H II region 
(Whiteoak & Gardner 1977) suggests that the uniform density assumption is not valid and 
that the bubble may be expanding asymmetrically into a region of lower density on that 
side. The derived expansion timescale should be considered a lower limit. 



5. Conclusions 

We have used Hubble Space Telescope WFPC2 imaging to obtain photometry for 1136 
stars in the Wd2 field, many of which are blends of multiple sources that are difficult to 
distinguish using ground-based observations. With Magellan and SOAR optical spectroscopy 
we identified fifteen new 0-type stars, including two new probable binary systems: #664 
(03 + 05.5) (MSP 188) and #714 (03 + ?). The photometry and spectral type information 
allowed us to determine the reddening parameters Ry and Ay individually for the 15 new 
O stars and 14 additional O stars from the literature. The mean Ry of 3.77 is larger than 
the canonical Galactic value and leads to a larger inferred extinction and a smaller distance 
modulus compared to most previous studies. Using Padova stellar isochrones we determined 
a new spectroscopic and photometric distance to Westerlund 2 of {d) = 4.16 ±0.06 (random) 
+0.26 (systematic) kpc. In archival Spitzer GLIMPSE images we discovered a faint mid-IR 
PAH ring that surrounds #889 (MSP 18) which seems to lie at the center of a secondary 
cluster located ~55" to the north of the cluster core. 

The distance revision of Wd2 calls into question whether the Wolf-Rayet star WR 20a 
and the TeV 7-ray sources HESS J1023-575 and HESS J1023-5746 are indeed cluster mem- 
bers or even located at a similar distance. Our distance determination allows that the HESS 
objects could be cluster members, but they would be almost twice as luminous when placed 



at 4.16 kpc instead of 2.4 kpc (Ackermann et al. 2011). The distance of 4.16 kpc potentially 



leaves WR 20a in isolation and calls into question ideas about the formation of massive WR 



stars. If WR 20a is at a distance of ~8 kpc like Rauw et al. (2005, 2007) and Rauw et al. 



(2011) conclude, then, did it form in isolation, and its relative position to Wd2 is just coin- 
cidence? It could be argued that neither the WR star nor the 7-ray sources belong to Wd2. 
However, if we use R07's photometry for WR 20a in conjunction with our Ry determination, 
then this star would lie at a distance of ~4.44 kpc, well within the uncertainties of our new 
distance, thus making its association with Wd2 likely. Our new distance determination also 
motivates a revision of Wd2's mass. A distance of 4.16 kpc will revise its mass downward 
compared to estimates using larger distances and will make Wd2 less comparable to the 
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more massive clusters of the Milky Way. 
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Table 1. Observations 



Position 


Number 


Filter 


Exposure Time 


UT 




of exposures 




(sec) 


Date 


1 


3 


F336W 


160 


2005-05-19 


1 


3 


F439W 


40 


2005-05-19 


1 


3 


F555W 


8 


2005-05-19 


1 


3 


F814W 


20 


2005-05-19 


2 


3 


F336W 


50 


2005-05-20 


2 


3 


F439W 


12 


2005-05-20 


2 


3 


F555W 


2 


2005-05-20 


2 


3 


F814W 


2 


2005-05-20 



Table 2. Photometry 



ID MSP91 a (J2000 °) S (J2000 °) F336W (tf336W F439W crF439W F555W (tfsssw F814W o-f8iw U B V I S/N Source 



1 




155 


93865 


-57 


75719 


99 


999 





000 


99 


999 





000 


17 


621 





023 


16 


585 





009 


99 


999 


99 


999 


17 


597 


16 


547 


147 


5 




2 




155 


94239 


-57 


75466 


17 


757 





016 


17 


640 





014 


16 


687 





009 


15 


521 





004 


17 


780 


17 


567 


16 


661 


15 


478 


320 







3 




155 


94329 


-57 


75762 


99 


999 





000 


99 


999 





000 


20 


203 





075 


17 


438 





010 


99 


999 


99 


999 


20 


270 


17 


461 


109 


1 




66 




155 


96852 


-57 


75353 


99 


999 





000 


99 


999 





000 


99 


999 





000 


20 


894 





094 


99 


999 


99 


999 


99 


999 


99 


999 


11 


4 




67 


273b 


155 


96877 


-57 


76713 


18 


633 





018 


18 


534 





016 


17 


529 





013 


16 


454 





006 


18 


656 


18 


462 


17 


503 


16 


411 


231 


7 




68 




155 


96880 


-57 


75063 


99 


999 





000 


99 


999 





000 


20 


986 





137 


18 


910 





022 


99 


999 


99 


999 


20 


996 


18 


890 


49 


7 




69 




155 


96893 


-57 


74199 


99 


999 





000 


99 


999 





000 


99 


999 





000 


21 


443 





111 


99 


999 


99 


999 


99 


999 


99 


999 


10 


2 




70 




155 


96899 


-57 


75914 


99 


999 





000 


99 


999 





000 


21 


375 





176 


19 


061 





023 


99 


999 


99 


999 


21 


401 


19 


053 


47 


8 




71 




155 


96914 


-57 


75185 


99 


999 





000 


99 


999 





000 


99 


999 





000 


20 


238 





055 


99 


999 


99 


999 


99 


999 


99 


999 


19 


9 




72 


273a 


155 


96941 


-57 


76712 


15 


801 





004 


15 


118 





004 


99 


999 





000 


99 


999 





000 


15 


438 


15 


059 


99 


999 


99 


999 


429 


6 




1134 




156 


05588 


-57 


74563 


18 


930 





105 


18 


966 





127 


18 


005 





081 


16 


861 





077 


18 


991 


18 


895 


17 


979 


16 


819 


42 


2 


s 


1135 




156 


05922 


-57 


74585 


15 


547 





030 


15 


330 





022 


14 


785 





030 


14 


082 





024 


15 


548 


15 


298 


14 


761 


14 


045 


264 


7 


s 


1136 




156 


05973 


-57 


74560 


19 


166 





223 


99 


999 





000 


99 


999 





000 


18 


597 





202 


19 


362 


99 


999 


99 


999 


18 


560 


10 


1 


s 



00 
to 



Note. — 



Table [2]is published in its entirety in the electronic edition of this journal. A portion is shown here for guidance regarding its form and content. 
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Table 3. 


New Spectral Classification 








ID 


MSP91 


V 


B - V 


EW(A5410) 


EW(A5876) 


EW(5410)/EW(5876) 


T 


Spectral 


Comment 






(mag) 


(mag) 


(A) 


(A) 




(K) 


Type 






Stars with only absorption lines 


137 


165 


15.591 


1.660 


0.97 ± 0.09 


0.30 ± 0.07 


3.20 ± 0.79 


43000 ± 2000 


04 


S 


505^ 


196a 


16.094 


1.325 


0.33 ± 0.11 


0.37 ± 0.09 


0.89 ± 0.36 


34000 ± 2000 


08.5 


S 


528^ 


229a 


15.841 


1.515 


0.55 ± 0.10 


0.53 ± 0.09 


1.05 ± 0.26 


35000 ± 2000 


08 


s 


548 


151 


14.522 


1.409 


0.94 ± 0.04 


0.31 ± 0.06 


3.00 ± 0.58 


43000 ± 1000 


04 


s 


584 


157b 


15.442 


1.346 


0.87 ± 0.06 


0.79 ± 0.05 


1.10 ± 0.10 


35000 ± 700 


08 


M 


640 


233 


16.234 


1.312 


0.66 ± 0.12 


1.05 ± 0.11 


0.63 ± 0.13 


32000 ± 800 


09.5 


M 


664a 


188a 


13.464 


1.538 


0.12 ± 0.01 


0.02 ± 0.01 


5.70 ± 3.32 


48000 ± 4000 


03 


M 


664b 


188a 


13.464 


1.538 


0.05 ± 0.01 


0.02 ± 0.01 


2.04 ± 1.18 


40000 ± 4000 


05.5 


M 


704 


175a 


14.059 


1.354 


0.95 ± 0.06 


0.15 ± 0.05 


6.24 ± 1.92 


48000 ± 2000 


04 


M 


769 


219 


16.576 


1.445 


0.42 ± 0.12 


0.85 ± 0.05 


0.50 ± 0.14 


31000 ± 800 


09.5 


M 


771 


167b 


15.501 


1.433 


0.63 ± 0.06 


0.63 ± 0.06 


1.00 ± 0.13 


35000 ± 900 


08 


M 


857 


444a 


13.869 


1.410 


1.13 ± 0.04 


0.45 ± 0.05 


2.52 ± 0.29 


42000 ± 800 


04.5 


M 


879 


235 


16.645 


1.469 


0.59 ± 0.12 


0.77 ± 0.09 


0.77 ± 0.18 


33000 ± 1000 


09.5 


M 


896=^ 


183c 


16.246 


1.438 


0.66 ± 0.08 


0.80 ± 0.07 


0.83 ± 0.13 


34000 ± 900 


08.5 


M 


903 


183b 


15.789 


1.409 


0.77 ± 0.10 


0.85 ± 0.08 


0.90 ± 0.14 


34000 ± 1000 


08.5 


M 


1004 


32 


15.323 


1.325 


0.53 ± 0.06 


0.74 ± 0.05 


0.72 ± 0.09 


33000 ± 600 


09.5 


s 


Stars believed to be late O/early B 


260^'*' 


201 


16.920 


1.324 


< -0.15 


0.43 ± 0.12 








S 


547 


28 


17.383 


1.298 


< 0.19 


0.79 ± 0.14 








S 


549 


44a 


15.562 


1.300 


< 0.04 


0.70 ± 0.29 








S 


556 




17.124 


0.900 


< -0.21 


0.53 ± 0.11 








S 


620 


96a 


16.086 


1.279 


< -0.07 


0.18 ± 0.06 








S 


878 




16.615 


1.291 


< 0.85 


2.17 ± 0.94 








M 


Field Stars 


1028 




14.355 


0.373 


0.15 ± 0.08 


-0.31 ± 0.09 






A 


S 


Stars unable to be identified 


463^'*' 


125a 


18.982 


0.871 












S 


1012^ 


24 


15.601 


1.225 












S 


1026*' 


20a 


16.784 


1.289 


< 0.19 










S 


1031*'''= 




15.308 


1.257 


< 0.11 


-0.93 ± 0.07 








s 



Note. — Comment column lists the source of the spectroscopic data, M stands for the Magellan Telescope while the S stands for 
the SOAR Telescope. Negative values of EW indicate an emission line. 

*'This spectrum also exhibits Ha emission. 

''This spectrum is a blend of two or more stars. The labeled one is the dominant source on the slit. 
"^This source exhibits a Mg II absorption line. 



Table 4. Reddening Summary 



ID Spectral Type Ry (mag) Ay (mag) 



CCM89 F04 FM07 FM09 Average CCM89 F04 FM07 FM09 Average 



137 


04V 


3 


.76 


± 


0, 


.09 


3, 


.99 


it 


0.07 


3.84 


it 


0.07 


3 


.79 


it 


0.43 


3.85 


it 


0. 


,10 


7.55 


it 


0.09 


7 


.51 


it 


0. 


,21 


7.41 


it 


0, 


,22 


7.49 


it 


0.86 


7.49 


it 


0.06 


178 


04V 


3 


.74 


it 


0, 


.03 


3, 


,95 


it 


0.02 


3.93 


it 


0.03 


3 


.79 


it 


0.17 


3.85 


it 


0, 


,10 


6.40 


it 


0.03 


6 


.35 


it 


0, 


,05 


6.38 


it 


0, 


,07 


6.37 


it 


0.29 


6.37 


it 


0.02 


395 


07.5V 


3, 


.73 


± 


0, 


.03 


3, 


,96 


ih 


0.02 


3.77 


it 


0.03 


4. 


.12 


it 


0.33 


3.89 


± 


0, 


,18 


7.06 


it 


0.03 


7 


.03 


± 


0, 


,06 


6.92 


± 


0. 


,07 


7.19 


± 


0.91 


7.05 


it 


0.11 


505 


08.5V 


3 


.62 


± 


0, 


.06 


3, 


,86 


it 


0.05 


3.71 


it 


0.06 


4, 


.01 


it 


0.35 


3.80 


it 


0. 


,17 


6.34 


it 


0.06 


6 


.33 


it 


0. 


,12 


6.36 


it 


0, 


,14 


6.46 


it 


0.57 


6.38 


it 


0.06 


528 


08V 


3, 


.81 


± 


0. 


.06 


4. 


,04 


it 


0.05 


3.99 


it 


0.05 


3 


.87 


it 


0.32 


3.93 


it 


0. 


,11 


7.00 


it 


0.06 


6 


.97 


it 


0. 


,13 


6.97 


it 


0. 


,14 


6.99 


it 


0.58 


6.98 


it 


0.02 


548 


04V 


3 


.64 


± 


0, 


.05 


3. 


,90 


it 


0.04 


3.76 


it 


0.04 


3 


.73 


it 


0.25 


3.76 


it 


0. 


,10 


6.46 


it 


0.04 


6 


.46 


it 


0. 


,09 


6.48 


it 


0. 


,10 


6.43 


it 


0.43 


6.46 


it 


0.02 


549^ 


Bl.OV 


4, 


.01 


± 


0, 


.05 


4, 


,32 


it 


0.03 


4.01 


it 


0.04 


3 


.91 


it 


0.22 


4.06 


it 


0. 


,18 


6.12 


it 


0.04 


6 


.11 


it 


0. 


,07 


6.09 


it 


0, 


,08 


6.03 


it 


0.35 


6.09 


it 


0.04 


584 


08V 


3, 


.69 


± 


0, 


.02 


3, 


,94 


it 


0.01 


3.73 


it 


0.02 


3 


.84 


it 


0.13 


3.80 


it 


0, 


,11 


6.22 


it 


0.02 


6 


.22 


it 


0, 


,04 


6.19 


it 


0, 


,05 


6.29 


it 


0.21 


6.23 


it 


0.04 


597 


06.5V 


3, 


.70 


± 


0, 


.04 


3, 


,95 


it 


0.03 


3.76 


it 


0.03 


3 


.80 


it 


0.20 


3.80 


it 


0, 


,11 


6.41 


it 


0.03 


6 


.40 


it 


0, 


,07 


6.39 


it 


0, 


,08 


6.42 


it 


0.34 


6.40 


it 


0.01 


620 


BIV 


3, 


.70 


± 


0, 


.04 


3, 


,93 


it 


0.03 


3.82 


it 


0.04 


3 


.76 


it 


0.23 


3.80 


it 


0, 


,10 


5.76 


it 


0.03 


5 


.74 


it 


0, 


,07 


5.77 


it 


0, 


,08 


5.71 


it 


0.35 


5.75 


it 


0.03 


640 


09.5V 


3 


.64 


± 


0, 


.02 


3, 


,87 


it 


0.01 


3.73 


it 


0.02 


4, 


.12 


it 


0.15 


3.84 


it 


0. 


,21 


6.38 


it 


0.02 


6 


.35 


it 


0. 


,04 


6.37 


it 


0, 


,05 


6.55 


it 


0.25 


6.41 


it 


0.09 


704 


04V 


3.68 


± 


0. 


.14 


3.92 


it 


0.11 


3.76 


it 


0.12 


3 


.75 


it 


0.72 


3.78 


it 


0. 


,10 


6.27 


it 


0.12 


6 


.27 


it 


0. 


,27 


6.27 


it 


0. 


,29 


6.28 


it 


1.20 


6.27 


it 


0.01 


714 


03V 


3, 


.64 




0. 


.06 


3. 


,83 


it 


0.04 


3.73 


it 


0.05 


3 


.61 


it 


0.29 


3.70 


it 


0. 


,10 


6.08 


it 


0.05 


6 


.04 


it 


0. 


,10 


6.08 


it 


0. 


,11 


6.01 


it 


0.49 


6.05 


it 


0.04 


722 


06V 


3 


.59 


± 


0.02 


3, 


,79 


it 


0.01 


3.65 


it 


0.01 


3 


.91 


it 


0.22 


3.74 


it 


0. 


,14 


7.24 


it 


0.02 


7 


.19 


it 


0. 


,04 


7.23 


it 


0, 


,04 


7.32 


it 


0.71 


7.25 


it 


0.05 


738 


05.5V 


3, 


.62 


± 


0, 


.04 


3, 


,85 


it 


0.03 


3.73 


it 


0.04 


3 


.71 


it 


0.22 


3.73 


it 


0, 


,09 


6.03 


± 


0.04 


6 


.01 


it 


0, 


,07 


6.02 


it 


0, 


,08 


6.05 


it 


0.36 


6.03 


it 


0.02 


769 


09.5V 


3, 


.62 


± 


0, 


.02 


3, 


,84 


± 


0.01 


3.65 


it 


0.02 


3 


.94 


it 


0.14 


3.76 


it 


0, 


,15 


6.66 


± 


0.02 


6 


.63 


it 


0, 


,04 


6.63 


± 


0, 


,06 


6.76 


it 


0.29 


6.67 


it 


0.06 


771 


08V 


3, 


.72 


± 


0, 


.03 


3, 


,94 


± 


0.01 


3.80 


± 


0.03 


3 


.86 


± 


0.15 


3.83 


± 


0, 


,09 


6.64 


it 


0.03 


6 


.61 


it 


0, 


,05 


6.60 


± 


0, 


,06 


6.65 


± 


0.27 


6.62 


it 


0.03 


804 


06V 


3 


.58 


± 


0, 


.02 


3. 


,79 


it 


0.01 


3.71 


it 


0.01 


3, 


.66 


it 


0.09 


3.69 


it 


0. 


,09 


6.84 


it 


0.02 


6 


.81 


it 


0. 


,03 


6.91 


it 


0, 


,04 


6.82 


it 


0.17 


6.85 


it 


0.05 


826 


09.5V 


3 


.60 


± 


0. 


.03 


3. 


,77 


it 


0.02 


3.67 


it 


0.03 


3 


.89 


it 


0.17 


3.73 


it 


0. 


,12 


6.54 


it 


0.03 


6 


.45 


it 


0. 


,06 


6.49 


it 


0, 


,07 


6.55 


it 


0.32 


6.51 


it 


0.04 


843'' 


04.5V 


3, 


.53 


± 


0. 


.05 


3. 


,79 


it 


0.04 


3.64 


it 


0.04 


3 


.70 


it 


0.24 


3.67 


it 


0. 


,11 


6.32 


it 


0.05 


6 


.33 


it 


0. 


,10 


6.32 


it 


0, 


,11 


6.42 


± 


0.41 


6.35 


± 


0.05 


857 


04.5V 


3 


.47 


± 


0, 


.03 


3, 


,65 


it 


0.02 


3.63 


it 


0.03 


3 


.55 


it 


0.19 


3.58 


it 


0. 


,08 


6.14 


it 


0.03 


6 


.11 


it 


0. 


,07 


6.13 


it 


0, 


,08 


6.11 


it 


0.33 


6.12 


it 


0.02 


879 


09.5V 


3, 


.58 


± 


0, 


.03 


3, 


,79 


it 


0.02 


3.70 


± 


0.03 


4, 


.02 


it 


0.35 


3.77 


it 


0, 


,19 


6.91 


it 


0.03 


6 


.87 


it 


0, 


,06 


6.98 


± 


0, 


,07 


7.04 


± 


0.99 


6.95 


± 


0.07 


896 


08.5V 


3, 


.55 


± 


0, 


.03 


3, 


,76 


it 


0.01 


3.69 


it 


0.02 


4, 


.00 


it 


0.15 


3.75 


it 


0. 


,19 


6.60 


it 


0.03 


6 


.58 


it 


0. 


,05 


6.58 


it 


0, 


,06 


6.78 


it 


0.27 


6.64 


it 


0.10 


903 


08.5V 


3, 


.58 


± 


0, 


.03 


3, 


,79 


it 


0.02 


3.67 


it 


0.03 


3 


.97 


it 


0.17 


3.75 


it 


0. 


,17 


6.71 


it 


0.03 


6 


.68 


it 


0. 


,06 


6.72 


it 


0, 


,07 


6.85 


it 


0.31 


6.74 


it 


0.08 


903 


06.5V 


3 


.56 


± 


0, 


.05 


3. 


,79 


it 


0.04 


3.66 


it 


0.04 


3 


.95 


it 


0.27 


3.74 


it 


0. 


,17 


6.73 


it 


0.05 


6 


.72 


it 


0. 


,11 


6.78 


it 


0, 


,12 


6.87 


it 


0.49 


6.77 


it 


0.07 


906 


06.5V 


3 


.53 


± 


0, 


.05 


3, 


,74 


it 


0.04 


3.64 


it 


0.05 


3 


.99 


it 


0.28 


3.72 


it 


0. 


,19 


6.56 


it 


0.05 


6 


.54 


it 


0. 


,11 


6.51 


it 


0, 


,12 


6.74 


it 


0.49 


6.59 


it 


0.11 


913 


03.5V 


3, 


.56 


± 


0, 


.05 


3, 


,80 


it 


0.04 


3.66 


it 


0.04 


3 


.61 


di 


0.26 


3.65 


it 


0, 


,10 


6.41 


it 


0.05 


6 


.40 


it 


0, 


,10 


6.42 


it 


0, 


,11 


6.36 


it 


0.46 


6.40 


it 


0.03 


924 


08V 


3, 


.51 


± 


0, 


.03 


3, 


,74 


it 


0.02 


3.60 


it 


0.03 


3 


.88 


± 


0.17 


3.68 


it 


0, 


,16 


6.47 


it 


0.03 


6 


.47 


it 


0, 


,06 


6.40 


it 


0, 


,07 


6.59 


it 


0.30 


6.48 


it 


0.08 


1004'' 


09.5V 


3, 


.65 


± 


0, 


.02 


3, 


,83 


it 


0.01 


3.79 


it 


0.02 


3 


.62 


± 


0.13 


3.72 


it 


0, 


,10 


6.00 


it 


0.02 


5 


.92 


it 


0, 


,03 


5.94 


it 


0, 


,04 


5.83 


it 


0.21 


5.92 


it 


0.07 


1039 


04.5V 


3 


.44 


± 


0, 


.04 


3, 


,63 


it 


0.03 


3.47 


it 


0.04 


3, 


.69 


it 


0.31 


3.56 


it 


0. 


,12 


6.43 


it 


0.04 


6 


.39 


it 


0. 


,09 


6.42 


it 


0, 


,10 


6.42 


it 


0.74 


6.42 


it 


0.02 



Table 4 — Continued 



ID Spectral Type 




Ry (mag) 










Av (mag) 




CCM89 


F04 


FM07 


FM09 


Average 


CCM89 


F04 


FM07 


FM09 Average 


3.64 ± 0.11 


3.86 ± 0.13 


3.73 ± 0.11 


3.84 ± 0.15 


3.77 ± 0.09 


6.51 ± 0.38 


6.48 ± 0.37 


6.49 ± 0.37 


6.54 ± 0.42 6.51 ± 0.38 



00 



''The IR, photometry of these stars comes from the 2MASS PSC instead of |Ascenso et al. [ ( 12007] . 
'^This star lacks J photometry. 



ID Spectral Type 



CCM89 



Table 5. Reddening Table for the HST Fi 

%336W (™^g) 



Filters 



%439W ('"^g) 



F04 



FM07 



Average 



CCM89 



F04 



FM07 



Average 



^F555W 



(mag) 



CCM89 



F04 



FM07 



Average 



137 


O4.0V 


11 


36 


± 


0.03 


11 


17 


± 


0.17 


11.24 


± 


0.18 


11.35 


± 


0.03 


9.80 


± 


0.03 


9.53 


± 





17 


9.54 


± 


0.18 


9.79 


± 





03 


7.68 


± 





02 


7.56 


± 


0.16 


7.54 


± 


0.17 


7.68 


± 


0.02 


178 


O4.0V 


9 


65 


± 


0.03 


9 


51 


± 


0.05 


9.53 


± 


0.06 


9.61 


± 


0.02 


8.32 


± 


0.03 


8.10 


± 





05 


8.09 


± 


0.06 


8.25 


± 





02 


6.52 


± 





02 


6.41 


± 


0.04 


6.43 


± 


0.06 


6.49 


± 


0.02 


395 


07.5V 


10 


67 


± 


0.03 


10 


51 


± 


0.05 


10.45 


± 


0.06 


10.61 


± 


0.02 


9.18 


± 


0.03 


8.95 


± 





05 


8.94 


± 


0.06 


9.11 


± 





02 


7.18 


± 





02 


7.09 


± 


0.05 


7.05 




0.06 


7.16 




0.02 


505 


08.5V 


9 


70 


± 


0.03 


9 


55 


± 


0.11 


9.57 


± 


0.12 


9.69 


± 


0.03 


8.31 


± 


0.03 


8.10 







11 


8.12 


± 


0.12 


8.29 


± 





03 


6.46 


± 





02 


6.38 




0.10 


6.38 


± 


0.11 


6.45 




0.02 


528 


O8.0V 


10 


49 


± 


0.03 


10 


35 


± 


0.11 


10.36 


± 


0.12 


10.48 


± 


0.03 


9.07 


± 


0.03 


8.84 


± 





11 


8.82 


± 


0.12 


9.05 


± 





03 


7.13 


± 





02 


7.03 


± 


0.10 


7.03 


± 


0.11 


7.12 


± 


0.02 


548 


O4.0V 


9 


86 


± 


0.03 


9 


69 


± 


0.08 


9.74 


± 


0.09 


9.83 




0.02 


8.45 


± 


0.03 


8.23 


± 





08 


8.24 


± 


0.09 


8.42 


± 





02 


6.58 


± 





02 


6.50 


± 


0.08 


6.49 


± 


0.09 


6.57 


± 


0.02 


549* 


Bl.OV 


8 


96 


± 


0.03 


8 


84 


± 


0.07 


8.99 


± 


0.07 


8.95 


± 


0.02 


7.83 




0.03 


7.63 


± 





07 


7.67 


± 


0.07 


7.79 


± 





02 


6.23 


± 





02 


6.15 


± 


0.06 


6.11 


± 


0.07 


6.22 


± 


0.02 


584 


O8.0V 


9 


43 


± 


0.03 


9 


28 




0.03 


9.34 


± 


0.04 


9.37 


± 


0.02 


8.11 


± 


0.03 


7.89 


± 





03 


7.91 


± 


0.04 


8.00 







02 


6.33 


± 





02 


6.24 


± 


0.03 


6.23 


± 


0.04 


6.29 


± 


0.01 


597 


06.5V 


9 


71 


± 


0.03 


9 


54 


± 


0.06 


9.60 


± 


0.07 


9.67 


± 


0.02 


8.35 


± 


0.03 


8.12 


± 





06 


8.14 


± 


0.07 


8.29 


± 





02 


6.52 


± 





02 


6.43 


± 


0.06 


6.41 


± 


0.07 


6.51 


± 


0.02 


620 


Bl.OV 


8 


73 


± 


0.03 


8 


59 


± 


0.06 


8.62 


± 


0.08 


8.70 


± 


0.02 


7.51 


± 


0.03 


7.30 


± 





06 


7.30 


± 


0.08 


7.46 


± 





02 


5.87 


± 





02 


5.77 




0.06 


5.77 




0.07 


5.85 


± 


0.02 


640 


09.5V 


9 


74 




0.03 


9 


58 


± 


0.03 


9.61 


± 


0.05 


9.67 


± 


0.02 


8.35 


± 
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Averages 



9.92 ± 0.00 9.87 ± 0.01 9.85 ± 0.01 9.90 ± 0.58 8.50 ± 0.00 8.37 ± 0.01 8.35 ± 0.01 8.46 ± 0.49 6.63 ± 0.00 6.59 ± 0.01 6.57 ± 0.01 6.62 ± 0.38 



Table 5 — Continued 



ID Spectral Type %336W (™^g) "^F439W (™^S) %555W (™^g) %814W (^^s) 

CCM89 F04 FM07 Average CCM89 F04 FM07 Average CCM89 F04 FM07 Average CCM89 F04 FM07 Average 

^The IR photometry of these stars comes from the 2MASS PSC instead of |Ascenso et aL| | |2007] . 
''This star lacks J photometry. 



Table 6. Distance Calculation. 



ID MSP91 Spectral F555W (F439W - F555W) -^FSSSW (F439W - F555W)o -^pggg-y^f DM d 

Type (mag) (mag) (mag) (mag) (mag) (mag) (kpc) 
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*These spectral classification were determined by R07 and Rll. 



Fig. 1. — Spitzer three-color image of RCW49 and the Wd2 cluster with 4.5 /im in blue, 5.8 
fim in green, and 8.0 fim in red. The dashed and solid outlines depict the fields covered by 
the "long" (Position 1) and "short" (Position 2) HST WFPC2 pointing, respectively. 
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Fig. 2. — Typical photometric uncertainty as a function of instrumental magnitude for each 
of the HST filters. Asterisks denote the "long" exposures (red in the electronic edition) and 
dots denote the "short" exposures (blue in the electronic edition). 
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Fig. 3. — Difference between the "short" and "long" exposure magnitudes versus long expo- 
sure magnitude for stars detected on both sets of images for each of the HST filters. The 
mean differences are consistent with zero, within the noise. The rms values are, 0.2162, 
0.1567, 0.1577, and 0.0887 magnitudes, respectively, slightly larger at fainter magnitudes. 
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Fig. 4. — Position map of the 1136 observed stars in Wd2. 



Fig. 5. — RGB image composite of Wd2 using the F439W (blue), F555W (green), and 
F814W (red) HST filters. Note that the probable member stars near the cluster core appear 
similarly red, while probable foreground stars having less reddening appear white. 
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Fig. 6. — Same as figure [5] for a zoomed view of the cluster core. Key stars to this investi- 
gation are labeled according to the numeration system created for this work. 



Fig. 7. — Same as Figure [5| showing a close-up view of the cluster core. Note the bright 
star in the upper right (not labeled) that appears less red than other stars, marking it as a 
foreground object. 
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Fig. 8. — Photometric differences between MSP91 and this work. The HST photometry 
is fainter than the ground-based observations. The numeration follows the identification 
system of MSP91. 
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Fig. 9. — Photometric differences between R07 and this work. Once again the HST pho- 
tometry is fainter than the ground-based observations. The numeration follows the systems 
of MSP91 and R07. 
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Fig. 10. — CMD of Wd2 using the HST instrumental filters. Asterisks (red colored in the 
online edition) represent the stars that were previously observed by MSP91. Key stars are 
labeled according to the numeration system created for this work. 
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Fig. 11. — Normalized spectra obtained at both the Magellan and SOAR telescopes. Each 
spectrum is offset by a constant value for easier visualization. Under each spectra is the stars 
identification number assigned in this work along with the spectral type. Spectra appear 
in order of descending temperature/spectral type from top to bottom. The main stellar 
absorption lines are labeled at the top of the figure (along with the interstellar Na I D lines), 
and the strongest diffuse interstellar band (DIB) hues are identified at the bottom the figure. 
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Fig. 12. — Continuation of Figure 



11 



- 51 - 




5000 



5500 6000 

A (A) 



6500 



Fig. 13. — Similar to Figure 11 but for the field stars or stars with indeterminate spectral 



type. The main stellar absorption lines are labeled at the top of the figure (along with the 
interstellar Na I D lines), and the strongest diffuse interstellar band (DIB) lines are identified 
at the bottom the figure. 
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Fig. 14. — Spitzer three-color image of the diffuse mid-IR ring surrounding the stellar sub- 
group centered on 7^889 (MSP 18). The image displays 4.5 /im in blue, 5.8 /im in green, 
and 8.0 fim in red. MSP 18 has been masked out for easier visualization of the ring. 




Fig. 15. — Velocity plot of He I A5876 (top panel), and He II A5411 (bottom panel) lines of 
the probable binary system ^^664 showing that each line appears to have two components 
separated by 360 - 490 km s~^. 
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Fig. 16. — MIKE spectrum of #714 showing the very weak He I, and strong He II indicative 
of a hot star. Note the residual nebular [O III] emission. 
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Fig. 17.— MIKE spectrum of #738, classified as 05. Note the residual nebular [O III] 
emission. 
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Fig. 18. — Example reddened magnitude SEDs for stars #137, 528, 704, and 1039. The 
solid lines are the reddened magnitudes, the squares are the HST flight photometry and the 
circles are the Ascenso et al. (2007) IR photometry. The solid line was reddened with the 
CCM89 reddening law, the dashed line used F04, and the dotted hne used FM07 (these are 
colored in blue, cyan, and magenta, respectively, in the electronic edition). 



-57- 



-0.6 



■ reddened ZAMS 



_ _> = 3.77, = 6.51 
curved track 



-0.5 



= 5.2^ 



^ -0.4 

CO 

I -0.3 
m 

^ -0.2 



-0.1 







' * ^ ^ ^ 
A = 6.03 * ^ ^) 

* * 

*• • • I ^ 

A -ii^*" 



* 



•A^ = 7.54 



0.1 
0. 



1.2 1.4 1.6 

F439W - F555W 



1.8 



Fig. 19. — CCD showing a section of the reddened ZAMS fit to the cluster Wd2. The aster- 
isks (red in the electronic edition) denote the stars having spectral types discussed in this 
work, R07, and Rll. The ZAMS was reddenened by Ry = 3.77, Ay = 6.51 with correspond- 
ing color excess of E(F336W - F439W) = 1.47, and E(F439W - F555W) = 1.81. Dashed 
and dotted lines represent a section of the linear and curved reddening tracks, respectively, 
for early-type stars (colored in magenta and cyan in the color edition). Numerical labels 
mark extinctions. Ay, along the curved reddening track. 
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Fig. 20. — Dereddened CMD of Wd2 at the adopted distance modulus of 12.98 mag. Labels 
show the IDs and respective spectral types. Labels at the right of the figure indicate the 
spectral types according to the ZAMS isochrone. The curves represent Padova isochrones, 
with ages as indicated in the legend. 



